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The Onset of Anisotropic Transport of Two-Dimensional Electrons in High Landau 
Levels: An Isotropic-to-Nematic Liquid Crystal Phase Transition? 

K. B. Cooper 1 , M. P. Lilly 1 !, J. P. Eisenstcin 1 , L. N. Pfeiffer 2 , and K. W. West 2 
1 California Institute of Technology, Pasadena CA 91125 
2 Bell Laboratories, Lucent Technologies, Murray Hill, NJ 07974 

The recently discovered anisotropy of the longitudinal resistance of two-dimensional electrons near 
half filling of high Landau levels is found to persist to much higher temperatures T when a large 
in-plane magnetic field -By is applied. Under these conditions we find that the longitudinal resistivity 
scales quasi-linearly with B^/T. These observations support the notion that the onset of anisotropy 
at B^— does not reflect the spontaneous development of charge density modulations but may 
instead signal an isotropic-to-nematic liquid crystal phase transition. 
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Transport experiments jl]|2| on clean two-dimensional 
electron systems (2DES) at high perpendicular magnetic 
field have recently revealed a new class of collective 
phases at low temperatures. These new phases, which 
have been observed only when electrons occupy the ex- 
cited Landau levels, are quite different from the fractional 
quantized Hall liquids found almost exclusively in the 
lowest (N = 0) Landau level j| . The new phases exhibit 
strong electrical transport anisotropies near half filling 
and curious re-entrant insulating behavior in the flanks 
of the third and higher (N > 2) Landau level (LL). At 
the phenomenological level, these properties are remark- 
ably consistent with those expected from the "stripe" and 
"bubble" charge density waves (CDW) predicted on the 
basis of Hartree-Fock (HF) theory . These inhomoge- 
neous states are stabilized by an exchange energy-driven 
tendency to phase separate which is particularly effective 
in the excited LLs. Near half filling, the ground state is 
approximately a unidirectional CDW, i.e. a stripe phase. 
The anisotropy of the longitudinal resistance observed at 
this filling is presumed to result from very different elec- 
trical conductivities parallel and perpendicular to stripes 
whose orientation relative to the crystal axes of the host 
semiconductor is determined by some, as yet unknown, 
symmetry-breaking field. 

With the filling factor v defined as the ratio of the 2D 
electron density n s to the degeneracy eB/h of a single 
spin resolved LL, the simplest HF stripe state consists 
of parallel strips in which the filling factor v alternates 
between adjacent integers. For example, at average fill- 
ing v — 9/2, corresponding to half filling of the lower 
spin branch of the N=2 LL, the local filling factor al- 
ternates between v = 4 and v — 5. The period of 
this alternation is set by the competition between the 
direct and exchange Coulomb interactions and is esti- 
mated to be on the order of 100 nm in typical samples. 
These same interactions also determine the temperature 
at which the stripes form and HF estimates are in the 
few Kelvin range 0]. This contrasts with the experi- 
mental observation that resistive anisotropy only sets in 



below about 100 mK. Although disorder in the 2D sys- 
tem may account for some of this discrepancy, a more 
interesting idea is that the onset of anisotropy does not 
in fact signal the initial development of charge density 
modulation, but instead reflects the orientational order- 
ing of local regions having pre-existing stripe order. This 
idea, first suggested by Fradkin and Kivelson , emerges 
naturally from the view that quantum and thermal fluc- 
tuations destroy the long range translational order of the 
stripe phases and render them analogous to nematic liq- 
uid crystals. The onset of anisotropy is therefore viewed 
as an isotropic-to-nematic phase transition. Wexler and 
Dorsey have estimated the transition temperature to 
be on the order of 200 mK, not far from the experimental 
result. In this paper, we report temperature-dependent 
magneto-transport studies which address this issue. Our 
results reveal that the sharp thermal onset of resistive 
anisotropy is replaced by a heavily smeared transition 
when a strong in-plane magnetic field component B« is 
added to the existing perpendicular field. This obser- 
vation, which is reminiscent of the behavior of a ferro- 
magnet in a strong symmetry-breaking field, supports 
the isotropic-to-nematic liquid crystal transition idea. In 
addition, we report an intriguing scaling of the resistive 
anisotropy with T/Bu which may prove important in de- 
veloping a quantitative understanding of the putative ne- 
matic phase. 

The sample used for the present experiments is a con- 
ventional single- interface GaAs/AlGaAs heterojunction 
grown by molecular beam epitaxy (MBE) on an [001]- 
oriented GaAs substrate. A thin sheet of Si impurities is 
positioned in the Alo.24Gao.76As 80 nm from the inter- 
face with the GaAs. Charge transfer from these donors 
results in a 2DES whose density and low temperature 
mobility are n s = 1.48 x 10 11 cm -2 and fj, = 1.1 x 10 7 
cm 2 /Vs, respectively. These parameters are obtained af- 
ter brief low temperature illumination of the sample by a 
red light emitting diode. The sample itself is a 4x4 mm 
square cleaved from the parent wafer. The edges of this 
square are parallel to the [110] and [110] crystallographic 
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FIG. 1. a) Longitudinal resistance anisotropy around 
v = 9/2 at T — 25 mK. Solid trace: R xx \ average current flow 
along [110] crystallographic direction. Dashed trace: R yy \ av- 
erage current flow along [110]. b) Temperature dependence of 
resistances at v — 9/2. c) R xx and R yy at v = 9/2 at T = 25 
mK vs. in-plane magnetic field along [110] and [110]. 



directions. Eight diffused In ohmic contacts are posi- 
tioned at the corners and midpoints of the sides of the 
sample. Longitudinal resistance measurements are per- 
formed by driving a 20 nA, 13 Hz a.c. current between 
opposing midpoint contacts and detecting the resulting 
voltage between corner contacts. We denote by R xx and 
Ryy the resistances obtained for average current flow par- 
allel to the [110] and [110] directions, respectively. Ad- 
justable in-plane and perpendicular magnetic fields are 
obtained by in situ tilting of the sample relative to the 
field B of a single superconducting solenoid Q . 

Figure la summarizes the resistance anisotropy effect 
observed around filling factor v = 9/2. The solid and 
dashed curves represent R xx and R yy , respectively, at 
T = 25 mK in a purely perpendicular magnetic field (i.e. 
•B||= 0). The anisotropy is largest at half-filling, where 
Rxx/ Ryy reaches a value of approximately 60 in this sam- 
ple. As Fig. lb shows, raising the temperature causes R xx 
to fall rapidly and R yy to rise somewhat. Above about 
60 mK both resistances are saturated at values less than 
60 ft. This transport anisotropy is also present at filling 
factors v = 11/2, 13/2, 15/2, etc., albeit with decreasing 
strength. As reported previously [|l],|j, no anisotropy is 
observed in the N = 1 first excited Landau level jll| nor 
in the N = lowest LL. 

An in-plane magnetic field can alter the orientation 
of the resistive anisotropy in high Landau levels PJL0||, 
Figure lc summarizes this effect by displaying the de- 
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FIG. 2. Temperature dependencies of longitudinal resis- 
tance in hard and easy directions at v = 9/2. a) and b): 
Rhard vs. T for various in-plane magnetic fields By directed 
along [110] and [110], respectively, c) and d): Analogous re- 
sults for R e as y . Note magnified vertical scale. 



along [110], R xx quickly falls to a small value just as R yy 
quickly rises to a large value: The in-plane field inter- 
changes the "hard" and "easy" transport directions. On 
the other hand, with B\\ directed along [110], no such in- 
terchange occurs. These data demonstrate that large in- 
plane fields can overwhelm the native symmetry-breaking 
potential in the sample and force the hard resistance di- 
rection to be parallel to the in-plane field. Jungwirth et 
al. [[l2| have argued that this effect arises from the finite 
thickness of the 2DES. For samples such as the present 
one, their HF theory reveals an energetic advantage for 
the stripes to lie perpendicular to the in-plane magnetic 
field. By comparing the in-plane field required to inter- 
change the anisotropy axes with this theory, an estimate 
of the native anisotropy energy may be obtained: Typical 
values are around 1 mK per electron [^3| . 

Figures 2a and 2b display the temperature dependen- 
cies of the longitudinal resistance in the "hard" and 
"easy" directions, Rhard and R easy , at v = 9/2. The 
data in Fig. 2a were taken with the in-plane magnetic 
field directed along [110]. In this case, the in-plane field 
does not interchange the principal axes of the resistance 
anisotropy and thus Rhard = Rxx- For Fig. 2b, however, 
the in-plane field is along [110] and therefore does inter- 
change the hard and easy directions, provided B\\ > 0.5 
T. With the exception of the Bn = data, the traces in 
Fig. 2b correspond to R yy , which is Rhard at the in-plane 



pendences of R xx and R vv at v = 9/2 and T = 25 mK fields chosen. Figures 2c and 2d show the corresponding 
on In the right half of the panel B\\ is directed along results for R easy , the longitudinal resistance in the easy 
the [110] direction; in the left half it is along [110]. For B\\ direction, on a magnified scale. Note that R easy is small 



2 



compared to Rhard for all Bu at low temperatures. 

It is clear from Figs. 2a and 2b that a large Bu sys- 
tematically broadens out the variation of Rhard at low 
temperatures. At B\\ = the bulk of the change in Rhard 
occurs in a relatively narrow window below T w 60 mK. 
In contrast, at Bu = 9 T Rhard diminishes smoothly up to 
temperatures above 600 mK. This qualitative broadening 
effect due to the in-plane magnetic field occurs for both 
orientations of the field relative to the crystal directions. 
Interestingly, the effect is not obviously present in R ea sy > 
although close inspection does reveal some softening at 
low temperatures. 

The thermal broadening of the resistive anisotropy in 
high Landau levels induced by strong in-plane magnetic 
fields is reminiscent of the behavior of a mean-field fer- 
romagnet in the presence of an external magnetic field 
fl4j| . Very weak external fields merely break rotational 
symmetry and orient the magnetization of the system 
which otherwise develops spontaneously below the Curie 
temperature T c . In the presence of a strong external field, 
however, the ferromagnetic transition is broadened out in 
temperature and substantial magnetization is present at 
temperatures T > T c . For T » T c the magnetization is 
that of the free spins in the system. By analogy, our ob- 
servation that the high Landau level resistive anisotropy 
persists to high temperatures in the presence of large B\\ 
suggests that local stripe "moments" exist in the 2DES at 
such temperatures. Lacking a strong symmetry-breaking 
in-plane field, these moments are disordered at high tem- 
perature and transport is isotropic. This scenario is con- 
sistent with the nematic-to-isotropic liquid crystal tran- 
sition proposed by Fradkin and Kivelson || . 

To pursue this idea further we have extracted, in an 
approximate manner, the microscopic resistivities of the 
2DES from the measured macroscopic resistances. As 
emphasized by Simon JlB) , the resistances Rhard and 
Reasy in square samples exaggerate the actual anisotropy 
of the underlying resistivities phard and p eaS y Jljl • This 
stems from the fact that the current distribution in an 
anisotropic sample is quite different from that in an 
isotropic one. In particular, currents flowing in the easy 
direction are channeled in toward the axis connecting 
the source and drain contacts. As a result, the volt- 
ages present at the remote contacts used to determine 
Reasy can be extremely small and therefore particularly 
sensitive to local inhomogeneities and other defects in 
the sample. Conversely, current flow in the hard direc- 
tion is spread out more uniformly across the sample and 
the voltages used to determine Rhard are robust. As a 
consequence, we will focus our attention on phard- 

Figure 3a displays the temperature dependence of 
Phard at v — 9/2 for several Bu along [110]. These re- 
sistivities were computed from the measured resistances 
Rhard and Reasy under the assumption of Simon's Jig ] 
classical model of the current distribution in the sample. 
While this model has shortcomings, its success in rec- 
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FIG. 3. a) Temperature dependence of phard at v = 9/2 for 
various Bu along [110]. b) Temperature dependence of order 
parameter of the 2D XY model for various symmetry-breaking 
fields h/ J. Arrow denotes Kosterlitz-Thouless transition tem- 
perature for h — 0. 



onciling the differences between the resistive anisotropy 
in square and Hall bar samples suggests that it provides 
reasonable estimates of the resistivities Jig ]. In any case, 
Fig. 3a reveals that phard exhibits the same broaden- 
ing of its temperature dependence as does Rhard when a 
large in-plane magnetic field is applied. 

Fradkin, et al. |]|l^] and Wexler and Dorsey |jj have 
argued that the stripe system undergoes a finite tem- 
perature Kosterlitz-Thouless (KT) transition from an 
isotropic state possessing only local stripe order to a ne- 
matic phase whose director field possesses algebraically 
decaying correlations. Below the KT transition temper- 
ature, the system is singularly susceptible to external 
symmetry breaking fields and is thus readily oriented by 
some (as yet unexplained fl3|| ) weak symmetry break- 
ing field in the GaAs/AlGaAs sample. Fradkin, et al. 
have modeled the system as a classical 2D XY magnet, 
modified to accomodate a director field. The order pa- 
rameter for the system is taken to be M =<cos(26*) >, 
where 9 is the angle between the director and an arbi- 
trarily weak symmetry-breaking field h. Using Monte 
Carlo methods, they have been able to achieve good fits 
to the observed temperature dependence of the resistive 
anisotropy at v = 9/2 (in a different sample) using a sen- 
sible value of the exchange energy, J, and a weak sym- 
metry breaking field, h ~ 0.05 J. We have extended these 
calculations to larger h to simulate the effect of the large 
in-plane magnetic field used in the present experiments 
p9| . As Fig. 3b shows, the temperature dependence of 
M is systematically broadened out as h increases and 
bears a strong qualitative resemblance to the behavior of 
Phard (and Rhard) at large B\\. 

We turn now to an intriguing scaling property exhib- 
ited by the longitudinal resistance Rhard- Figures 4a and 
4b display the same Rhard data plotted against the scaled 
temperature T / Bu (for non-zero Bu). Plotted in this way, 
the data collapse onto a single curve with impressive pre- 
cision. Figs. 4c and 4d show log-log plots of the computed 
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FIG. 4. a) and b): Rhard vs. the scaled temperature T/Bu 
for B|| directed along [110] and [110], respectively, c) and d): 
Log-log plot of the hard axis resistivity, phard, calculated from 
the measured resistances, vs. T/B^. Dashed lines are power 
laws: phard'X- (T/B\\)~ a with a = 0.75. 



resistivities phard vs. T/Bu and again the scaling is evi- 
dent. This result demonstrates that the hard resisitivities 
depend on temperature and in-plane field only through 
the ratio T/B\\ . The log-log plots reveal that this depen- 
dence is quasi-linear: The dashed lines in the figure are 
proportional to (T/ B\\)~ - 75 . We note in passing that 



the easy direction resistance R eaS y does not exhibit this 
scaling behavior. Although the origin of this difference is 
not known, it may be related to the aforementioned cur- 
rent channeling effects which are inherent to transport in 
anisotropic conductors. 

The connection between the scaling behavior of Rhard 
and the nematic-to-isotropic phase transition picture is 
not yet understood. While at sufficiently high temper- 
ature T solutions to the 2D XY model must eventually 
approach a Curie law M ~ h/T, this limiting behav- 
ior is not well developed in the {h, T) parameter range 
which seems most relevant to our experimental data. On 
the other hand, a quantitative comparison of the 2D XY 
model with our data is hindered by several factors. For 
example, the relation between the symmetry-breaking 
parameter h and the in-plane magnetic field Bu is sub- 
tle and highly sample-specific [l2| . Also, Bu may affect 
the transport mechanism of current along and across the 
stripes to some degree. Even the appropriateness of the 



classical 2D XY model itself can be questioned. We em- 
phasize, however, that the experimental observation of 
scaling is robust and offers a new insight into the nature 
of the anisotropic phases. 

In the picture outlined above, the essential role of the 
in-plane magnetic field is to orient previously existing lo- 
cal striped regions. It is, however, also possible that Bu so 
alters the energetics of the electron system that it creates 
stripe order at temperatures far above where anisotropic 
transport appears spontaneously. If this is the case, then 
the nematic-to-isotropic transition interpretation of our 
data may be unjustified. Although we cannot completely 
rule out such a scenario, there are good reasons to dis- 
count it. First, the data shown above clearly indicate 
that the in-plane field severely broadens the onset of re- 
sistive anisotropy; it does not merely shift the sharp tran- 
sition encountered at Bu = to a higher temperature. 
Second, the scaling behavior we observe is sub-linear in 
T/Bu. A model in which the both stripe formation and 
orientation are dependent on Bu would presumably lead 
to a super-linear dependence on Bu . 

In conclusion, we have found the anisotropy character- 
istic of transport in half-filled high Landau levels to be 
heavily broadened by a strong in-plane magnetic field. 
Although we have concentrated here on the v = 9/2 
state, the same effect is observed at v = 11/2 and other 
half- filled high Landau levels. This finding suggests that 
local stripe moments exist at relatively high tempera- 
tures where the transport is isotropic in the absence of 
the in-plane field. In addition, we have also observed a 
remarkable scaling of the resisitive anisotropy with the 
ratio T/B\\. 
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